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Abstract 



^Sj ' A constraint on masses of superheavy gauge and Higgs multiplets at the grand 

^ ' unification (GUT) scale is obtained from the gauge couphng unification in the case of 

high-scale supersymmetry. We found that all of the particles may lie around a scale 
^T) \ of 10^^ GeV so that the threshold corrections to the gauge coupling constants at 

^^ ■ the GUT scale are smaller than those in the case of the low-energy supersymmetry. 

■^ . In addition, the GUT scale tends to be slightly lower when the gauginos are heavier 

and, thus, the proton decay rate via the A-boson exchange process is expected to 

be enhanced. 



1 Introduction 

Supersymmetric grand unified theories (SUSY GUTs) [T],|2] are promising candidates of 
pliysics beyond tlie Standard Model (SM). Indeed, tfiey are strongly motivated by an 
experimental observation which implies that the gauge coupling constants of the SM gauge 
groups are to be unified at a certain high-energy scale with good accuracy [3H7] . From the 
observation, the unified scale is estimated as Mgut — 2 x 10^^ GeV. In the GUTs, new 
superheavy particles are supposed to appear around the scale, which make the couplings 
run together above the scale. These particles are naturally expected to have masses of 
the order of Mgut- They are, of course, beyond the reach of collider experiments, and 
there is little hope to search for them or to measure their masses directly. 

In Refs. [SlIS] , a way of constraining the masses of superheavy particles indirectly by re- 
quiring the gauge coupling unification is discussed and limits on the masses are presented. 
Later, by applying the same method with more accurate gauge coupling constants, the 
authors in Ref. [lOj derive a more stringent constraint on the masses of the color-triplet 
Higgs boson, the adjoint Higgs bosons, and the X bosons in the context of the minimal 
SUSY SU(5) GUT [HIS]- They have found that while the masses of the adjoint Higgs and 
the X bosons are around 10^^ GeV, the mass of the color-triplet Higgs boson lies in the 
region of 3.5 x 10^^ GeV < Mhc < 3.6 x 10^^ GeV, which is significantly below the GUT 
scale. This implies that the threshold corrections to the gauge coupling constants are still 
not small, even in the SUSY GUTs. In fact, the analysis is quite sensitive to the mass 
spectrum in the intermediate scale, such as those of the SUSY particles. In Ref. [TU], 
all of the particles except for gauginos are assumed to be at 1 TeV. The gaugino masses 
are set to be around the electroweak scale with the GUT relation for the gaugino masses 
being assumed. 

Currently, on the other hand, the SUSY models with heavy sfermions have been widely 
discussed with a lot of attention [TTHTB]. Although such models originally have been dis- 
cussed as a possible candidate of SUSY models from a theoretical point of view, now they 
are also supported by the latest results of the LHC experiments; no significant excess in 
the SUSY searches |T9H21] and the discovery of the 126 GeV Higgs boson [22||23] both sug- 
gest that the SUSY breaking scale is somewhat higher than the electroweak scale. While 
this high-scale SUSY scenario is difficult to be probed at the collider experiments, it still 
has a lot of phenomenological consequences which might be checked in other experiments. 
Recent studies on the subject are given in Refs. |23H36]. 

As mentioned above, the GUT scale mass spectrum inferred from the indirect analysis 
presented in Refs. |8l|9] is highly dependent on the SUSY mass spectrum. In particular, if 
the SUSY particles have masses much larger than 0{1) TeV, previous results are expected 
to be changed significantly. In this Letter, therefore, we revisit the analysis in the high- 
scale SUSY scenario. We carry out the calculation in the minimal SUSY SU(5) GUT 
with sfermions having masses well above the electroweak scale. We will see that while 
the constraint on the adjoint Higgs and X boson masses only differs from previous ones 
slightly, that on the color-triplet Higgs mass is found to be changed by more than an order 
of magnitude and actually improved in the sense that it may also be around the GUT 



scale. Interestingly enough, this result again indicates that the high-scale SUSY scenario 
is rather supported than the traditional low-energy SUSY models. 

This Letter is organized as follows: in Sec. |2l the high-scale SUSY model which we 
discuss below is presented, and its phenomenological aspects are briefly described. The 
mass spectrum of the superheavy particles in the minimal SUSY SU(5) GUT is also 
displayed there. In the subsequent section, we discuss the method of constraining the 
masses of the GUT-scale multiplets by means of the renormalization group equations 
(RGEs) as well as the threshold corrections of the gauge couplings. Then, we show some 
results in Sec. HI Section. [5] is devoted to conclusions and discussion. 

2 Model and Spectrum 

Let us begin by presenting a high-scale SUSY model discussed in this Letter. We consider 
the particle content of the minimal supersymmetric Standard Model (MSSM). Then, the 
only assumption which we adopt here is that there exists a djTiamical SUSY-breaking 
sector with the Kahler potential having a generic structure. Then, all of the scalar bosons 
except the lightest Higgs boson acquire masses of the order of the gravitino mass 7/13/2, 
while the lightest Higgs boson mass is flne-tuned to be rrih ~ 126 GeV. The higgsinos may 
in general have similar masses to the gravitino mass, though they might be suppressed 
if there are some extra chiral symmetries. The gaugino masses are, on the other hand, 
generated by the quantum effects pTl[37] and, thus, suppressed by a loop factor compared 
with m3/2. There are two kinds of contributions to the gaugino masses; one is the anomaly 
mediation effect pT|IB7] and the other is from the higgsino-Higgs boson loop diagram. 
These two effects give rise to the gaugino masses as 

33 ^? / 1 

M2 = -^(m3/2 + L) , 

Ms = -3^^3/2 , (1) 

with 

L = /i^ sin 2/3 , f^ ^ In ^ (2) 

representing the higgsino-Higgs boson loop contribution. Here, Ma and ga (a = 1,2, and 
3) are the U(l)y, SU(2)2,, and SU(3)c gaugino masses and gauge coupling constants, 
respectively. We use the SU(5) normalization for the U(l)y coupling, i.e., gi = ^/5J3g'. 
Further, fi^ and rriA denote the higgsino and the heavy Higgs boson masses, respectively, 
which are assumed to be the same order of magnitude as sfermion masses, tan /3 is the 
ratio of the vacuum expectation values (VEVs) of the Higgs flelds; tan (3 = (Hu) / {Ha) . 
Since the higgsino mass is presumed to be around the gravitino mass, the higgsino-loop 
contribution L is also expected to be as large as 7713/2- The values of the gaugino masses 



are, however, dependent on the relative phase between fi^ and m3/2. In the following 
discussion we just assume the gaugino masses are lighter than the scalar masses by a loop 
factor, and regard them as free parameters. 

As mentioned in Introduction, this model has a lot of fascinating features from a 
phenomenological point of view. They are originally discussed in Refs. [TIHT8] , and recent 
development is given in Refs. [23H3S] • In these works, a typical scale for sfermion masses 
is taken to be C(IO^-IO^) TeV, which explains the 126 GeV Higgs boson mass. With such 
heavy particles the SUSY flavor and CP problems [38] as well as the gravitino problem 
are considerably relaxed. In this case, the gaugino masses are 0{1) TeV because of an 
one-loop factor. Note that from Eq. ([1]) it is found that with a moderate value of L, wino 
becomes the lightest among gauginos, and thus the lightest SUSY particle in this model. 
It is quite interesting since the wino dark matter with a mass of 2.7-3 TeV is consistent 
with cosmological observations [39] . The dark matter might be searched directly [28ll40p41] 
or indirectly [26lll2lll3] in future dark matter experiments. 

Moreover, the gauge coupling unification in this high-scale SUSY model is found to 
be achieved as precisely as that in the MSSM. Thus the SUSY GUT is still promising in 
the case of high-scale supersymmetry. In the next section, we in turn require the gauge 
coupling unification and constrain the GUT scale mass spectrum by using the requirement. 

Now, we briefly summarize the superheavy particles in the minimal SUSY SU(5) 
GUT, which we adopt as the working hypothesis in this paper. The SUSY SU(5) gauge 
theory includes twenty-four gauge superfields V"^ with A = 1, ... ,24. By using the SU(5) 
generators T^ we define a 5 x 5 matrix representation of the vector superfields such that 
V = V^T^, with the components written as 
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We collectively call Xa and Ya the X-bosons hereafter. The unified gauge group SU(5) 
is spontaneously broken by the VEV of the adjoint Higgs boson S^ {A = 1, . . . , 24) to 
SU(3)c;xSU(2)ixU(l)y. Again, we define E = S^T^ and its components by 
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The MSSM Higgs superfields are, on the other hand, incorporated into fundamental and 
anti-fundamental representations as follows: 
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where 
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(6) 



are the MSSM Higgs superfields. Hq and Hca are called the color-triplet Higgs multiplets. 
The superpotential of the Higgs sector in the minimal SU(5) SUSY GUT is given by 



W = ^AsTrE^ + ^msTrS^ 



XhH^H + muHH . 



(7) 



After the adjoint Higgs field gets the VEV, (S) = 1/ ■ diag(2, 2, 2, -3, -3) {V = ms/As), 
the SU(5) gauge group is broken to the SM SU(3)cxSU(2)ixU(l)y gauge groups without 
breaking the supersymmetry. Also the parameter rriH is fine-tuned as mn = 3XhV in 
order to realize the doublet-triplet mass splitting in H and H. With the condition, the 
mass of the color-triplet Higgs boson is Mhc = Mjj^ = 5XhV. The X-boson mass is 
Mx = 5\/2g5V with g^ the unified gauge coupling constant. As regards the adjoint Higgs 
multiplets, the components S3 and Eg have masses of M^ = M^g = M^g = |AsV^ and E24 
has Ms24 = |AeI^. The components S(3* 2) and S(3 2) become the longitudinal component 
of the X-bosons, and thus do not show up as physical states. 



3 Renormalization Group Analysis 

In this section, we present the RGEs of the gauge and Yukawa coupling constants, as well 
as the boundary conditions at each threshold. We use the DR scheme [H] in this work. 
First, we write down the two-loop beta functions [IS]. In the MSSM, the two-loop RGEs 
for the gauge coupling constants are given as 
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(10) 



with i/i {i = t, b, t) the top, bottom and tau Yukawa coupling constants, respectively. Since 
the Yukawa couplings enter into the two-loop level contributions to the gauge coupling 
RGEs, it is sufficient to consider the RGEs for the Yukawa couphngs at one-loop level. 



They are given as 



d 
d 
d 



IGvr^ 

1 
IGvr^ 

1 
16^ 



Vt 



Vb 



Gy^ 


+ yl- 


LO 2 


-^gl 


Gyl 


+ y^t+yl- 


7 2 


^yl 


+ ^yl 


9 2 


-H 



16 ' 

- 3^3 


5 


9 16 9 

3^2 - y^3 



(11) 



Below the SUSY breaking scale {Ms)-, the squarks and sleptons, the higgsinos, and the 
heavy Higgs boson masses are decoupled so that the theory is regarded as the SM with 
gauginos. The contribution of gauginos and the SM particles to the coefficients of the 
beta functions is given as 
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(14) 



where the subscripts "SM" and "gaugino" indicate that the contributions are of the SM 
particles and gauginos, respectively. The running of the Yukawa couplings in this case is 
given as follows: 
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Next, we consider the matching conditions at each threshold scale. At the GUT scale, 
the gauge coupling constants in the SU(3)cxSU(2)iXU(l)y gauge theories are equated 
to the unified coupling constant g^ with the following threshold corrections at one-loop 
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Here, the conditions do not include constant (scale independent) terms since we use the 
DR scheme for the renormalization |5l[48] . From the equations it immediately follows that 
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The relations allow us to evaluate the masses of the heavy particles, Mh^ and MxMs, 
from the gauge coupling constants determined in the low-energy experiments through the 
RGEs [319]. While the couplings are well measured with high precision, the estimation is 
quite dependent on the spectrum in the intermediate region, especially on the masses of 
gauginos and higgsinos. 

For the SUSY breaking threshold, we just equate the gauge couplings above and below 
the threshold, and change the beta functions appropriately for each region. This approx- 
imation is valid since the particles appearing at the scale are assumed to be degenerate 
in mass. In the case of gauginos, on the other hand, we need to consider the threshold 
corrections since the mass difference among gauginos might be sizable. The condition is 
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(18) 



where 5'a(At)sM are the couplings in the SM while 5'a(Ai)gaugino are those above the gaugino 
threshold. 

The Yukawa couplings are matched as usual, i.e., at the SUSY breaking scale, the 
Yukawa couplings yi{fi) below the SUSY breaking scale are matched with the supersym- 



metric ones, yi{^)MSSM, as follows: 

yt{Ms)MSSM = ^—ryt{Ms) , 
snip 

yb{Ms)MSSM = pytiMs) , 

cosp 

yr{Ms)MSSM = oVAMs) . (19) 

Before concluding this section, we solve the RGEs at one-loop level and, taking the 
threshold corrections into account, derive relations between the superheavy masses and 
the low-energy gauge coupling constants. Such relations reflect the dependence of Mh^ 
and MxMy^ on the mass spectrum of the SUSY particles. By inserting to Eq. ( IT71) the 
one-loop solutions of the RGEs for the gauge couplings, we have 
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From Eq. (!20|) we flnd that the mass of the color-triplet Higgs Mh^ gets larger as the SUSY 
breaking scale Ms is taken to be higher. This originates from the mass difference among 
the components of the fundamental Higgs multiplets, i.e., the triplet-Higgs, higgsinos, 
heavy Higgs bosons, and the lightest Higgs boson. Therefore, the behavior of M^^ with 
respect to the SUSY breaking scale is universal in a sense. Further, Mh^ depends only on 
the ratio of M2 and M3. MxM^, is, on the other hand, independent of the SUSY braking 
scale Ms while dependent on the scale of the gauginos, not their ratio. This is because the 
right-hand side of Eq. f l2T]) results from the mass difference in the gauge vector multiplets 
and the adjoint Higgs multiplet, a part of which is included as the longitudinal mode of 
the gauge multiplets. It is also found that Mj^M^. decreases when the gaugino masses are 
taken to be large values. This is owing to the opposite sign of the contribution of gauge 
flelds to the gauge beta functions to those of matter flelds. This feature is, therefore, again 
model- independent. In the subsequent section, we carry out a similar analysis using the 
two-loop RGEs. 

4 Results 

Now we present some results for the RGE analysis which we discuss in the previous section. 
As noted above, the running of the gauge couplings is computed at two-loop level and the 
threshold corrections are taken into account at one-loop level. The masses of sfermions, 
heavy Higgs bosons, and higgsinos are taken to be Ms for brevity. Gaugino masses are 
assumed to be lighter than Ms by an one-loop factor. 

First, we consider the color-triplet Higgs mass Mhq- In Fig. [H we plot the dependence 
of Mhc on the SUSY breaking scale Ms in the pink lines. Here, the wino mass M2 is flxed 
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Figure 1: Predicted color-triplet Higgs mass Mh^ as functions of the SUSY breaking scale 
Ms (pink lines). Here, wino mass M2 is fixed to be 3 TeV and tan/3 = 3. Gluino-wino 
mass ratio, M3/M2, is set to be M3/M2 = 3, 9, and 30 from top to bottom, respectively. 
Theoretical errors coming from the strong coupling constant asimz) = 0.1184(7) ^9] 
are also shown. Horizontal blue line shows a result in the case of low-energy SUSY 
{Ms = 1 TeV, M2 = 200 GeV, and M3/M2 = 3.5). 
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Figure 2: Color-triplet Higgs mass Mh^;, as functions of gluino mass M3 (pink lines). Here, 
tan/3 = 3 and Ms = 10^ TeV. Upper and lower lines correspond to M2 = 3 TeV and 
300 GeV, respectively. Error bars indicate the input error of the strong coupling constant 
o^si'iTT'z) = 0.1184(7) |1H]- Horizontal blue line shows a result in the case of low-energy 
SUSY {Ms = 1 TeV, M2 = 200 GeV, and M3/M2 = 3.5). 
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Figure 3: GUT scale Mqut = {M'^My^Y^^ as functions of gluino mass M3 (pink lines). 
Here, tan/3 = 3 and Ms = 10^ TeV. Upper and lower lines correspond to M2 = 300 GeV 
and 3 TeV, respectively. Error bars indicate the input error of the strong coupling constant 
Q^si'iTT'z) = 0.1184(7) [49j. Horizontal blue line shows a result in the case of low-energy 
SUSY (Ms = 1 TeV, M2 = 200 GeV, and M3/M2 = 3.5). 



to be 3 TeV, which is favored from the thermal relic abundance, and tan/3 = 3. The ratio 
of the gluino and wino masses, M3/M2, is set to be M3/M2 = 3, 9, and 30 from top to 
bottom, respectively. Further, we show the error of the calculation coming from that of 
the strong coupling constant as{mz) = 0.1184(7) [Hj. The horizontal blue line shows a 
result in the case of low-energy SUSY {Ms = 1 TeV, M2 = 200 GeV, and M3/M2 = 3.5) 
as a reference. In this case, we have 8.6 x 10^^ < Mh^ ^ 1-4 x 10^^ GeV. This figure 
well illustrates the feature read from the approximated expression given in Eq. ( 12U]) : M^c 
increases as the SUSY breaking scale grows while it decreases when the ratio M3/M2 
becomes large. To see the latter feature more clearly, we show its dependence on the 
gluino mass M3. Again, we set tan/3 = 3, and the SUSY breaking mass is fixed to be 
Ms = 10'^ TeV. The upper and lower lines correspond to M2 = 3 TeV and 300 GeV, 
respectively. These two figures show that M^c is strongly dependent on Ms and M3/M2. 
Therefore, to predict the mass with high accuracy, precise determination of the masses 
of gauginos as well as the SUSY breaking scale is inevitable. Any way, in the high-scale 
SUSY scenario it is found to be possible for the mass of the color-triplet Mhc to be around 
~ 2 X 10^^ GeV, which is expected by the gauge coupling unification. 

Next, we discuss constraints on MjrMj] derived from the relation (1211) . From now 
on we define Mqut = (M'^My^Y^^ and refer to it as the GUT scale. The equation (1^ 
tells us that the GUT scale depends on only the gaugino masses at the leading order, so 
we express Mgut as functions of the gaugino masses. In Fig. [3] we plot it as functions 
of gluino mass. Here again we fix tan/3 = 3 and Ms = 10^ TeV. The upper and lower 



lines correspond to M2 = 300 GeV and 3 TeV, respectively. Again, the horizontal blue 
line shows a result in the case of low-energy SUSY with Ms = 1 TeV, M2 = 200 GeV, 
and M3/M2 = 3.5, which gives Mqut — 1-9 x 10^^ GeV. The error bars indicate the 
input error of the strong coupling constant as{mz) = 0.1184(7) [39], though the effect 
is negligible. We see that the GUT scale also has little dependence on the gaugino 
masses. In that sense, the prediction is robust compared with that for M^c- However, 
as discussed in the previous section, the GUT scale Mqut gets lower when the gaugino 
masses become larger (Mqut oc (M3M2)~^'^). This feature is quite interesting when one 
considers the proton decay via the X-boson exchange processes. Although the change in 
Mgut is small, it might be significant since the proton decay lifetime scales as oc Mjr. 
For instance, when Mx — 0.8 x 10^^ GeV, the proton lifetime via the X-boson exchange 
reduces to around 5 x 10^^ years jSU], which is slightly above the current experimental 
bound, t{p -^ e+7r°) > 1.29 x 10^^ yrs [51]. 

Finally, we briefly comment on the tan/3 dependence of the results. Although the 
one-loop computation is not related with tan/3, the two-loop results might be affected 
through the running of the Yukawa couplings. We have found, however, that the effects 
on the results are negligible. 

5 Conclusions and Discussion 

In this Letter, we have presented constraints on the masses of the GUT scale particles 
from the gauge coupling unification in the case of the high-scale SUSY scenario. To that 
end, we have used the two-loop RGEs for the gauge couplings with one-loop threshold 
corrections considered. As a result, the mass of the color-triplet Higgs multiplets M^c 
turns out to be considerably large compared with previous results in the traditional low- 
energy SUSY scenario, while the GUT scale is found to be slightly lower. These are 
generic features resulting from the mass difference among the components of the same 
supermultiplet of the SU(5) gauge group. Interestingly, all of the superheavy particles 
might be around 10^^ GeV in the high-scale SUSY models. 

The mass spectrum of the GUT scale particles predicted here stimulates us to recon- 
sider the proton decay in the case of high-scale supersymmetry. As mentioned to above, 
the relatively low GUT scale enhances the proton decay rate via the X-boson exchang- 
ing process. If the enhancement is strong enough, the proton decay in the p — )■ e"'"7r° 
channel might be searched in future experiments. Furthermore, the experiments may also 
reach the proton decay through the color-triplet Higgs exchange. Since the decay process 
predicts too short lifetime [10], some suppression mechanism for the process has been 
assumed. Nevertheless, with Mhc larger than those considered in the previous litera- 
ture and sfermions much heavier than the electroweak scale, this process might evade the 
current experimental bound without any mechanism of limiting the color-triplet Higgs 
exchanging process. In such a case, it is possible for the p — )■ K~^iJ mode, which is the 
main decay mode in the case of the color-triplet Higgs exchange, to be searched in near 
future. A detailed analysis of this decay process will be given elsewhere 
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Note that, in proposed models where the color-triplet Higgs exchange is suppressed by 
some mechanism, large threshold corrections to the gauge coupling constants at the GUT 
scale tend to appear. One of the examples is introduction of the Peccei-Quinn symmetry 
[53] . It was found that when the threshold corrections at the GUT scale are small, 
the suppression mechanism does not work [Ml [55] . In the SUSY SU(5) GUTs in higher 
dimensional space the U(l)ij symmetry forbids the dimension- five proton decay, though 
the Kaluza-Klein particles generate large threshold corrections to the gauge coupling 
constants [56j. The high-scale supersymmetry would be another solution for the proton 
decay problem, in which such large threshold corrections are not required. 

Acknowledgments 

The work of N.N. is supported by Research Fellowships of the Japan Society for the 
Promotion of Science for Young Scientists. The work of J.H. is supported by Grant-in- 
Aid for Scientific research from the Ministry of Education, Science, Sports, and Culture 
(MEXT), Japan, No. 20244037, No. 20540252, No. 22244021 and No. 23104011, and 
also by World Premier International Research Center Initiative (WPI Initiative), MEXT, 
Japan. 

References 

[1] S. Dimopoulos and H. Georgi, Nucl. Phys. B 193, 150 (1981). 

[2] N. Sakai, Z. Phys. C 11, 153 (1981). 

[3] S. Dimopoulos, S. Raby and F. Wilczek, Phys. Rev. D 24, 1681 (1981). 

[4] W. Marciano and G. Senjanovic, Phys. Rev. D 25, 3092 (1982). 

[5] M.B. Einhorn and D.R. Jones, Nucl. Phys. B 196, 475 (1982). 

[6] U. Amaldi, W. de Boer and H. Furstenau, Phys. Lett. B 260, 447 (1991). 

[7] P. Langacker and M. -x. Luo, Phys. Rev. D 44, 817 (1991). 

[8] J. Hisano, H. Murayama and T. Yanagida, Phys. Rev. Lett. 69, 1014 (1992). 

[9] J. Hisano, H. Murayama and T. Yanagida, Nucl. Phys. B 402, 46 (1993) . 
[10] H. Murayama and A. Pierce, Phys. Rev. D 65, 055009 (2002) . 
[11] G. F. Giudice, M. A. Luty, H. Murayama and R. Rattazzi, JHEP 9812, 027 (1998) . 
[12] 
[13 



J. D. Wells, hep-ph/0306127 



N. Arkani-Hamed and S. Dimopoulos, JHEP 0506, 073 (2005). 



11 



[14] G. F. Giudice and A. Romanino, Nucl. Phys. B 699, 65 (2004) [Erratum-ibid. B 706, 
65 (2005)]. 

[15] N. Arkani-Hamed, S. Dimopoulos, G. F. Giudice and A. Romanino, Nucl. Phys. B 
709, 3 (2005). 

[16] J. D. Wells, Phys. Rev. D 71, 015013 (2005) . 

[17] L. J. Hall and Y. Nomura, JHEP 1003, 076 (2010) . 

[18] L. J. Hall and Y. Nomura, JHEP 1201, 082 (2012) . 

[19] G. Aad et al. [ATLAS Collaboration] , ■ar2(ivil20M949| [hep-ex] . 

[20] [ATLAS Collaboration], ATLAS-CONF-2012-109. 

[21] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. Lett. 109, 171803 (2012). 

[22] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 1 (2012). 

[23] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716, 30 (2012). 

[24] G. F. Giudice and A. Strumia, Nucl. Phys. B 858, 63 (2012). 

[25] M. Ibe and T. T. Yanagida, Phys. Lett. B 709, 374 (2012). 

[26] M. Ibe, S. Matsumoto and T. T. Yanagida, Phys. Rev. D 85, 095011 (2012). 

[27] L. E. Ibanez and I. Valenzuela. larXiv:1301. 51671 [hep-ph]. 

[28] J. Hisano, K. Ishiwata and N. Nagata, Phys. Rev. D 87, 035020 (2013). 

[29] K. S. Jeong, M. Shimosuka and M. Yamaguchi, JHEP 1209, 050 (2012). 

[30] R. Saito and S. Shirai, Phys. Lett. B 713, 237 (2012). 

[31] R. Sato, S. Shirai and K. Tobioka, JHEP 1211, 041 (2012). 

[32] B. Bhattacherjee, B. Feldstein, M. Ibe, S. Matsumoto and T. T. Yanagida, Phys. 
Rev. D 87, 015028 (2013). 

[33] L. J. Hall, Y. Nomura and S. Shirai, JHEP 1301, 036 (2013). 

[34] N. Arkani-Hamed, A. Gupta, D. E. Kaplan, N. Weiner and T. Zorawski, 
arXiv:1212.6971 [hep-ph]. 

[35] T. Moroi and M. Naeai. [arXiv:1303.0668l [hep-ph]. 

[36] D. McKeen, M. Pospelov and A. Ritz. arXiv: 1303. 11721 [hep-ph]. 

[37] L. Randall and R. Sundrum, Nucl. Phys. B 557, 79 (1999). 



12 



[38] F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvestrini, Nucl. Phys. B 477, 321 
(1996) . 

[39] J. Hisano, S. Matsumoto, M. Nagai, O. Saito and M. Senami, Phys. Lett. B 646, 34 
(2007). 

[40] J. Hisano, K. Ishiwata and N. Nagata, Phys. Lett. B 690, 311 (2010) . 

[41] J. Hisano, K. Ishiwata and N. Nagata, Phys. Rev. D 82, 115007 (2010). 

[42] J. Hisano, S. Matsumoto and M. M. Nojiri, Phys. Rev. Lett. 92, 031303 (2004). 

[43] J. Hisano, S. Matsumoto, M. M. Nojiri and O. Saito, Phys. Rev. D 71, 063528 (2005). 

[44] W. Siegel, Phys. Lett. B 84, 193 (1979). 

[45] M. E. Machacek and M. T. Vaughn, Nucl. Phys. B 222, 83 (1983). 

[46] S. Weinberg, Phys. Lett. B 91, 51 (1980). 

[47] L. J. Hall, Nucl. Phys. B 178, 75 (1981). 

[48] I. Antoniadis, C. Kounnas and K. Tamvakis, Phys. Lett. B 119, 377 (1982). 

[49] J. Beringer et al. [Particle Data Group Collaboration], Phys. Rev. D 86, 010001 
(2012). 

[50] J. Hisano, D. Kobayashi and N. Nagata, Phys. Lett. B 716, 406 (2012). 

[51] H. Nishino, K. Abe, Y. Hayato, T. lida, M. Ikeda, J. Kameda, Y. Koshio and 
M. Miura et al. ' arXiv: 12 03.40301 [hep-ex]. 

[52] J. Hisano, D. Kobayashi, T. Kuwahara, and N. Nagata, in preparation. 

[53] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977). 

[54] J. Hisano, H. Murayama and T. Yanagida, Phys. Lett. B 291, 263 (1992). 

[55] J. Hisano, T. Moroi, K. Tobe and T. Yanagida, Phys. Lett. B 342, 138 (1995). 

[56] L. J. Hall and Y. Nomura, Phys. Rev. D 64, 055003 (2001). 



13 



